Phase diagram of a d-wave superconductor with Anderson impurities 
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We present a self-consistent solution for a model of a d-wave superconductor with finite concen- 
tration of Anderson impurities at zero temperature using the slave boson method. We show how 
the phase diagram depends on the strength of interaction between impurity and extended states. 
For fixed impurity level energy Eo in the Kondo limit there is one superconducting-normal state 
transition for all impurity concentrations n. When Eo is close the Fermi energy there are three such 
transitions for impurity concentration exceeding certain minimum value n(To,Eo). If hybridiza- 
tion To is fixed and the depth of the impurity level is varied, there are always two transitions for 
concentration above n(Fo, Eo). 



The problem of magnetic and nonmagnetic impurities 
in superconductors (SC) has a long history. Recent years 
brought new experimental results, most of which concern 
strongly correlated systems such as high-temperature 
and heavy-fermion superconductors X Many of these 
studies are complicated by the uncertainty about the 
symmetry of the order parameter in various compounds 
as well as the precise nature of the interaction between 
localized and extended states. The proximity to antifcr- 
romagnetism in many of these materials poses additional 
challenges. 

Scanning tunneling microscopy experiments in 
Bi2Sr2CaCu208+,5 reveal inhomogeneities on a 
nanoscopic scale.— ~— Local variations of chemical 
composition, e.g. presence of excess oxygen atoms and 
cationic disorder outside the CuC>2 planes might be 
responsible for these modulations of the local gap and 
local density of states (DOS). In bismuth strontium 
copper oxide (BSCCO) doped with magnetic impurity 
Ni resonances in local density of states attributed to 
impurity were observed only in small-gap domains and 
no such resonances were seen in large-gap regions. 5 
This sensitivity of local DOS to the size of the gap 
may result from the proximity to the critical point 
where magnetic impurities become decoupled from the 
superconductor^— 

Motivated by these works we will try to understand 
how the phase diagram of a superconductor with 
magnetic impurities depends on the position of the 
impurity level, the strength of interaction between lo- 
calized state and conduction electrons and the impurity 
concentration. By changing the depth of the impurity 
level we can switch from the Kondo limit to the mixed 
valence regime. The phase diagram in these two cases 
is qualitatively different. Knowing the general form of 
the phase diagram may improve our understanding of 
experiments in impurity-doped superconductors. 



I. MODEL 

The model consists of electrons in a conduction 
band with BCS-type pairing interacting with Anderson 
impurity, 
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The operator c\ m creates electron in a spin-orbit partial 
wave state of angular momentum m and momentum k. 
The energies €k lie in a band of half-width D and con- 
stant density of states No = 1/2D. The impurity state 
has energy Eq and its hybridization matrix element with 
extended states is V. The constraint nj + b^b = 1 is 
added to prevent double occupancy of the impurity site. 

We assume a two-dimensional d-wave order parameter 
of the form A(fc) = Aocos(2(/)), where 4> is the angle in 
the k x — k y plane. In the mean field approximation the 
dynamics of the boson fields is neglected and the boson 
fields are replaced by their expectation values < >=< 
b >= z 1 / 2 . Minimizing the free energy with respect to 
the auxiliary boson field and position of the many-body 
resonance e/ we arrive at the mean field equations 
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where the conduction electron Green's function in clean 
superconductor G° is given by 

(G H)" 1 = W To-e fe T3-A(fc)(T 1 - l T 2 ) , (4) 
and the full impurity Green's function is 

Gj 1 ^) — lotq - e/r 3 - £/(w) . (5) 

Conduction electron Green's function is averaged over 
impurity positions in the usual way. We have to solve 
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FIG. 1: The phase diagram of d-wave superconductor with 
Anderson impurities for four values of impurity concentration. 
The lines separate the normal state from the superconductor 
(SC). The broken line shows location of the impurity quantum 
phase transition. Below that line the impurities are decoupled 
from the superconductor. The slope of this line is approxi- 
mately -0.26 which agrees with the numerical renormalization 
group^ result and a single-impurity large- N calculation for a 
d-wave superconductor 



the system of equations ©, together with the gap 
equation 
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For each frequency uj we solve self-consistently the 
Dyson equations for conduction electron and impurity 
self-energies 



E(w) = G°(w) 
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where n is the concentration of impurities, N is the de- 
generacy of the impurity energy level Eq, T = zTq = 
zirNoV 2 . Here we assumed N = 2. 

The self-energy in eqn. (|5|) in general con- 

tains both diagonal and off-diagonal terms in particle- 
hole space. Here the off-diagonal term is zero due to 
vanishing of the Fermi surface average of the off-diagonal 
part of the conduction electron Green's function. 

The renormalized frequencies are solutions of the fol- 
lowing equations 
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FIG. 2: The amplitude of the order parameter as a function 
of bare hybridization To for Eq — —0.1 and several impurity 
concentrations. The approximate location of the critical point 
where impurities are decoupled from the conduction band is 
marked with a circle. 



and 
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where brackets denote average over the Fermi surface. 
The presence of full Green's functions under integrals of 
eqs. (3), (4) and (7) means eqs. (10) and (11) have to be 
solved at each step in the integration routines. 



II. RESULTS 

The calculations were performed for nondegerate im- 
purities, N = 2, with superconducting order parameter 
A = 0.01. All energies are given in units of D. 

The obtained phase diagram is shown in Fig. 1. In the 
Kondo limit the superconducting-normal state boundary 
is a straight line. In the mixed valence regime the bound- 
ary for fixed Eq is nonmonotonic. 

Fig. 2 shows Ao as a function of To for Eq = —0.1. 
There is a local gap minimum at intermediate values of 

r , where Tk ~ T c0 , where Tk is defined as ^Jv 2 + ej . 

Pair breaking is strongest when the two energy scales 
are comparable. For n > 0.0117 there appears a normal 
state around Tq c ~ 0.04. For Tq < Tq c , in the region 
marked "A" in Fig. 2 the order parameter quickly rises 
to the value of the clean limit. Here the system is very 
sensitive to details of interaction. The steep curve at 
the boundary of region A means that small change of Tq 
or impurity concentration may completely alter system 
properties. The normal state is separated from the clean- 
like superconductor by only about 10% change of To. 
Experiments conducted in this regime may be subject to 
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FIG. 3: Critical concentration for the superconducting- 
normal transition as a function of To = -rNqV 2 for Eo = — 1 
and Eo = —0.1. The Eo = —1 curve shows the behavior in 
the Kondo limit. 



greater errors than that in other parts of the phase dia- 
gram. Also the presence of local disorder around each im- 
purity significantly influences the system behavior in this 
limit. Inhomogeneities in experimental samples, whether 
intrinsic to the superconducting state or influenced by 
impurity doping, may lead to situation where the system 
is normal in some regions and superconducting in the rest 
of the sample. 

For r > r 0c the impurity resonance widens and pair 
breaking at low energies is less effective. Further increase 
of Tq raises pair-breaking for all energies and eventually 
the normal state becomes more stable. Numerical diffi- 
culties prevent us from finding the exact location of the 
point where A / A c i oa n — > 1 • However this is not nec- 
essary for the purpose of this paper. We can see in Figs. 
3 and 4 that decoupling of impurities from the conduc- 
tion band occurs at the same critical value of r / E for 
different choices of Eo- 

The phase diagram in the n vs. r /|-Eo| plane is shown 
in Fig. 2. In the well developed Kondo limit, when the 
bare impurity energy level Eq is deep below the Fermi 
surface and impurity occupation number nj ~ 1 , there is 
only one superconducting-normal transition. This is il- 
lustrated by the Eg = — 1 curve in Fig. 3. The transition 
line is monotonic in the entire phase diagram. 

When Eo lies closer to the Fermi surface, see curve 
Eq = —0.1, there exists a critical concentration n' c such 
that for n < n' c the system remains superconducting for 
small and moderate values of ro/|-Eo| and the transition 
to the normal state occurs at r /|i?o| 3> 1 (beyond the 
range shown in Fig. 3). However for n > n' c there are two 
additional transitions for intermediate values of To/ 1 Sol, 
when Tk ~ T c q. 

This is clearly shown in Fig. 4 where Eo = —0.01. The 
horizontal axis of Fig. 4 is extended to larger values of To- 
This is the general form of the phase diagram when Eq 
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FIG. 4: The phase diagram in the mixed valence limit, Eo = 
—0.01. The scale of the horizontal axis is changed relative to 
Fig. 2 to illustrate behavior at large To- 
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FIG. 5: Superconducting order parameter Ao as a function 
of Eo for several values of To at n = 0.01. The unmarked 
dash-dotted curve is the result for r = 0.08. The vertical 
dotted line marks the impurity quantum phase transition. 

is close to the Fermi energy. The superconducting state 
boundary for Eq = —0.1 has similar shape (not shown 
in Fig. 3 due to smaller horizontal scale). For increas- 
ing Eo the phase transition line is shifted towards higher 
impurity concentrations and the normal state section at 
intermediate To/lSol becomes narrower. This can be ex- 
plained by the weaker pair breaking in the mixed valence 
regime. 

The impurity-induced peak in the conduction electron 
DOS remains at uj = for all Eo < 0. The peak splits in 
two, one at positive ui and one at negative lu when e/ 3> 
r. For a nondegenerate impurity in a superconductor 
with order parameter having lines of nodes this occurs 
only at So > 0. This conclusion was verified numerically. 
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The dependence of the order parameter on Eq for sev- 
eral values of r is shown in Fig. [5] This behavior is 
always nonmonotonic. The minimum Ao occurs when 
Tk ~ Ao- For larger To we have two superconducting 
regions, one in the Kondo limit and one in the mixed 
valence regime, separated by the normal state. 

III. CONCLUSIONS 

The impurity quantum phase transition at finite cou- 
pling is associated with the particle-hole asymmetry of 
the models If the impurity degeneracy N is increased, 
the impurity transition occurs at lower coupling. 12 This 
is expected since for larger N the impurity resonance 
is located at higher energy. The precise location of the 
main features of the phase diagram changes but its over- 
all qualitative form remains intactji 2 - 

The suppression of the superconducting state along the 
diagonal of the phase diagram in Fig. 1 results from the 
interplay of energy scales and is not directly related to 
the symmetry of the order parameter. In an s-wave su- 
perconductor the steepest initial decrease of the super- 
conducting critical temperature as a function of impurity 
concentration, (dT c / ' dn)\T c =T a0 occurs for T c q ~ Tjc*& At 
T = this corresponds to Ao ~ Tk and we expect a nor- 
mal state insertion similar to that in Fig. 1 also for order 
parameters without nodes on the Fermi surface. 

Features of the phase diagram should be visible in 
experiments by (a) changing impurity concentration at 
fixed carrier doping and (b) varying the strength of cou- 
pling to the impurity, the ratio Tq/Eq, at fixed impurity 
concentration. 

The depth of the impurity level may be tuned by apply- 
ing pressure. There exist studies of heavy-fermion com- 
pounds, most notably CeCu2(Sii_ x Ge x )2 where varying 
hydrostatic pressure reveals two superconducting regions 



in the phase diagram*^. The change of pressure shifts the 
chemical potential. Position of the bare impurity level Eq 
of Ce ions relative to Ef changes accordingly. This may 
explain the existence of two superconducting regions in 
the phase diagram of CeCu2(Sii_ x Ge x )2 similarly to re- 
sults of this work. Increasing pressure brings the system 
into the mixed- valence regime. It is equivalent to increas- 
ing the ratio To/I-Eol in our work. The existing explana- 
tion for the superconducting region at high pressure uses 
the concept of the valence-fluctuation mediated pairing 
mechanism^ It is interesting to note that a qualitatively 
similar nonmonotonic behavior of superconducting criti- 
cal temperature is obtained in an impurity problem when 
increasing Eq. We hope that similar experiments under 
varying pressure may be conducted on superconductors 
with Anderson impurities. They should reveal the details 
of the phase diagram described in this article. 

A nontrivial question is realistic treatment of finite im- 
purity concentrations in this problem. In real samples 
there are local inhomogeneities and variations of impurity 
concentration. There is some evidence that inhomogene- 
ity of the superconducting state may be intrinsic in some 
high-T c compounds even in absence of impurities in CuC>2 
planes. 16 This may have important consequences near the 
transition line in the region A of the phase diagram in 
Fig. 2, where a small change in hybridization or impu- 
rity level Eq implies a large change in properties. In this 
limit the entire sample might even consist of a mixture 
of normal and superconducting patches. This sensitivity 
should be also visible in the low-energy physics near the 
impurity critical point. In inhomogeneous samples pock- 
ets of clean-like superconductor may coexist with regions 
where low-energy physics is dominated by Kondo impu- 
rities. Scanning tunneling microscopy (STM) measure- 
ments should be helpful in investigating the phase dia- 
gram in the vicinity of the superconducting state bound- 
ary and the impurity critical point. 
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